Deformation and scattering in graphene over substrate steps 
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The electrical properties of graphene depend sensitively on the substrate. For example, recent 
measurements of epitaxial graphene on SiC show resistance arising from steps on the substrate. 
Here we calculate the deformation of graphene at substrate steps, and the resulting electrical 
resistance, over a wide range of step heights. The elastic deformations contribute only a very 
small resistance at the step. However, for graphene on SiC(OOOl) there is strong substrate-induced 
doping, and this is substantially reduced on the lower side of the step where graphene pulls away 
from the substrate. The resulting resistance explains the experimental measurements. 



The advance of very high speed graphene electronic^^ 
depends on understanding and controlling the interaction 
of graphene with the supporting substrate. Electron mo- 
bility can vary over many orders of magnitude depending 
on the substrate^^. Among other factors, morphologi- 
cal deformations of the graphene may limit mobility™^. 
It is therefore important to determine how the substrate 
morphology affects transport in supported graphene. 

Epitaxial graphene on SiC provides an ideal system 
in which to study the role of substrate morphology. 
SiC is a promising substrate because, in contrast to 
other approaches, it allows growth of epitaxial graphene 
directly on an insulating substrate^. However, epi- 
taxial graphene on SiC substrates generally exhibits 
smaller carrier mob ilities than exfoliated graphene on 
Si0 2 substratesPE^. The reason for this difference 
is not fully understood, but SiC substrates exhibit a 
high density of multilayer steps, which are implicated 
in the lower mobility. Several experiments show that re- 
sistance increases with step density^, step heights^ and 
step bunching}^, and the local electrical resistance asso- 
ciated with individual substrate steps has recently been 
measurecP, by scanning potentiometry in a scanning tun- 
neling microscope. 

Here we study graphene over an abrupt substrate step, 
as illustrated in Fig.[T] calculating both the structural 
deformation and the resulting electrical resistance. The 
results are directly applicable to epitaxial graphene on 
SiC, and also show more generally how the morphology 
affects electrical transport. We find that very little re- 
sistance arises directly from the structural deformations, 
despite the strong curvature of graphene as it passes over 
a step. For SiC, we nevertheless find a substantial resis- 
tance associated with the step, in good agreement with 
experiment. This resistance arises almost entirely from 
the electrical coupling between the graphene and sub- 
strate, which varies sharply in the vicinity of the step. 
Thus morphology plays a qualitatively different and far 
more important role in substrates such as SiC that dope 
the substrate or otherwise couple strongly, than it does 
for substrates such as Si0 2 that are electrically passive. 

We begin by determining the graphene geometry as 
it passes over a substrate step. The graphene defor- 
mation is determined by a balance between the van der 




FIG. 1: Illustration of graphene over a substrate step. Here 
h s is the step height and h eq is the equilibrium distance of 
graphene from the substrate surface due to van der Waals 
interaction, id is the length of graphene detachment from the 
substrate. 



Waals interaction, which favors conforming to the sub- 
strate, and elasticity, which favors keeping the graphene 
flat and smooth. Since the displacement field can vary 
on the atomic scale, we use an atomistic valence force 
model (VFM) to describe the elastic deformationifS]. 
The van der Waals interaction between graphene and 
substrate is modeled with the Lennard- Jones (LJ) 6- 
12 potentials. The parameters for our LJ model are 
determined by setting the equilibrium distance between 
graphene and substrate to h eq s» 3.4 A 16 , and the bind- 
ing energy to Eb ~ 40 meV per atoirP2E0. The total en- 
ergy is then simply the sum of these two contributions, 
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We calculate the minimum-energy geometry for 
graphene over a wide range of step height h s , allowing 
the graphene to slide to relax any in-plane strain. The 
relaxed geometries h(x) for three step heights typical of 
SiC are plotted in Fig. [2^. The presence of the atomic 
step increases the graphene area compared to when the 
substrate is flat. The extra length increases with h s and 
was found to be 0.1, 0.3, and 0.6 nm respectively. The 
nonlinear dependence reflects the increasing distortion 
(steeper slope) with increasing step height. The maxi- 
mum slope is of order 1, confirming the need for a fully 
numerical treatment. 
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FIG. 2: (a) Graphene geometry h(x) for various step heights as 
indicated. Symbols are the full numerical calculation. Solid lines 
are best fit using Eq.^ The respective step profile are illustrated by 
dashed lines, (b) Corresponding curvature k for geometries shown 
in (a), with lines and symbols as in (a). The minimum radius of 
curvature (inversion of maximum k) is indicated. 



Wc find that the calculated geometries can be well ap- 
proximated by a simple error function, 
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As shown in Fig.[2]D, even the variation in curvature 
across the step is well described by this simple functional 
form. The only noticeable discrepancy is that the error 
function is symmetric, while in the full calculation there 
is a slight asymmetry, with the maximum curvature in- 
duced at the upper edge of the step. For step height 
h s = 1.5 nm, we find a maximum curvature equivalent to 
that of a carbon nanotube of diameter 1.5 nm. 

The maximum curvature K max as a function of step 
height h s is plotted in Fig. [3^,, for both the lower and 
upper edge of the step. It is proportional to h s in the 
small h s limit. For the lower edge, the limit of large 
h s corresponds to the well-known problem of peeling 18 . 
In this limit, n max approaches y // 2Es /a/3 w 1.2 nm -1 , 
where a is area per carbon atom and fi s» 2.1 eV^is the 
bending rigidity. 

Figure [3]d summarizes the dependence of the graphene 
deformation h(x) on step height h s , in terms of the pa- 
rameters of Eq.[T] For a given step height, the charac- 
teristic step width d s determines the maximum radius of 
curvature r = 1/K, max - Another relevant lengthscale is 
the length id over which the graphene is detached from 
the substrate. For concreteness we define id to be the 
length of graphene separated by 2h eq or more from the 
substrate surface. Fig.[3]3 shows that remains zero at 
small h s but begins to increase rapidly for h s > h eq . At 
larger step height we find that id ~ 1.2 h s , which proves 
important in later discussions. 

Geometry -induced scattering: Curvature results in 
bandstructure changes that can scatter electrons near the 
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FIG. 3: (a) Numerically calculated maximum curvature Kmax for 
the upper (solid) and lower (open) edge of the step as function of 
step height h s . Note that the often used approximation for curva- 
ture, ft RJ d^h, is only reasonable for h s of few angstrom. Over the 
range of h s of our interests, it grossly overestimates the real cur- 
vature due to the large gradients in our graphene geometries, (b) 
Dependence of graphene deformation on step height h B . Graphene 
step width d s and lateral shift x a are obtained by fitting Eq.^to 
the relaxed geometry h(x) obtained numerically. We also show the 
detachment length i.e. the length of graphene separated from 
the substrate by 2h eq . Note that t& K, l.2h s for large h s . The step 
width d s depends only weakly on step height. 



step. To examine this effect, after performing the geome- 
try relaxation, we construct the Hamiltonian % within a 
nearest-neighbor Slater-Koster parameterized sp 3 tight- 
binding modeP^. We then calculate the transmission 
and the electrical resistanc e 1Z u sing the non-equilibrium 
Green's function formalis m 1 20 ! 21 ! in the limit of small volt- 
age across the step and no inel astic s cattering at the step. 
We use the known Fermi leve l 16 1 22 1 of Ef = 0.45 eV for 
graphene on SiC (0001). 

Since the maximum curvature increases with step 
height h s , the resistance also increases. For a step height 
of 1.5 nm, we obtain a resistance ~ 0.01 tt-nm. Figure[4^, 
includes results of a recent experiments^, which employed 
scanning potentiometry in a scanning tunneling micro- 
scope to resolve the resistance in graphene due to indi- 
vidual substrate steps. The measured resistance lZ exp has 
roughly linear dependence with step height h s , ~ 10 f2- 
/im for each nanometer step height. Evidently, the re- 
sistance due to curvature alone cannot account for this 
large lZ exp . While ir-a hybridization can result in new 
scattering states in the vicinity of Dirac poin1>22l, this ef- 
fect is significant only when r < 5 A221. Even for a very 
high step of h s = 1.5 nm, we find that the minimum ra- 
dius of curvature only shrinks to r as 5 A if we assume a 
much stronger van der Waals attraction, with a binding 
energy 80 meV, which seems unphysical. 

In our calculations thus far, we have ignored possi- 
ble in-plane strain inhomogeneities, which is known to 
result in electron scatterin g ^ 24 * 25 !. Due to the different 
thermal expansion coefficient between graphene and SiC, 
graphene can acquires a residual biaxial strain upon cool- 
ing if sliding is suppressed. Then due to the nonplanar 
geometry, graphene at the step could experience a uni- 
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FIG. 4: (a) Resistance through graphene due to substrate step, TZ, 
as function of step height h s . Various electrostatic models consid- 
ering doping variations are plotted, as described by Eq. [5] (Model 1) 
and Eq.|3l (Model 2). Experiment dataSis plotted in symbols, with 
statistical error bar indicated, (b) Potential profile Vd(x) derived 
from Eq.[2] (Model 1) for h = 0.5, 1.0, 1.5 nm. Red dashed line 
indicates the Fermi level position, (c) Same as (b), but derived 
from Eq. [3] (Model 2). The small hump observed for Vd(x) near the 
top is an artifact of the polynomial form of A c (h). 



axial stress relative to the rest of the sheet. Graphene on 
SiC substrates is reported to have strains from 0.1 — 1% 
as measured by Raman spectroscopy^. To estimate a 
very conservative upper bound for 72 due to strain inho- 
mogeneity, we consider a local tensile strain of 1% along 
the detached region, with the step along the zigzag di- 
rection where the scattering effect is largest!^!. The re- 
sult is 72, < 1 fi-/im. Thus some source of scattering 
much stronger than the geometrical deformations must 
be present to account for the measured resistance. 

Electrostatic doping effects: It is well known that when 
contacting graphene with metals, a difference in work- 
function results in electrostatic doping^H2U_ i n the case 
of SiC (0001), a similar doping occurs from the carbon 
buffer layer, w hich h as a high density of weakly dispersive 
interface states^^. This can be described by a capacitor 
model including quantum capacitance^, 

n(x) = 7 ( W sg h(x)n{x) — hvf^Trn(x) ] (2) 

where n{x) is the electron density in graphene, W sg is 
the workfunction difference between the carbon buffer 
layer and graphene, h is the distance between them, and 
7 is the buffer-layer density of states. We de note this 
as "Model 1" . In view of the flat bandfPS we take 



the limit of large 7 and we adjust W sg to reproduced 
the kno wn do ping at h = h eq , corresponding to heavy 
n-doping p 6 1 22 1 , with a Fermi level Ej — 0.45 eV. The ver- 
tical displacement h(x) changes the capacitive coupling, 
leading to doping variations. Substituting the relaxed 
geometry h(x) into Eq. 2j we calculate these variations. 
The associated potential shifts Vd{x) are shown in Fig.[4jD 
for different step heights. Increasing step height leads to 
larger detachment and doping variations. This translates 
to an increased 72. as shown in Fig.|4^,, still somewhat 
smaller than reported experimentally, but far larger than 
the scattering mechanisms previously discussed. 

In studies of metal ind uced doping of graphene, 
Khomyakov and co-worker d 28 * 29 ! reported that Eq.[2] 
could not properly describe the ab initio calculations of 
graphene on metals, presumably due to quantum me- 
chanical effects such as wavefunction overlap and corre- 
lations They suggested that t he ac curacy of Eq.[5]can be 
improved by the modificatio n 1 28 ! 29 ! 



h(x) -t h(x) + h* 

W sg + A c (h) 



W, 



(3) 



with h* and A c (h) approximated as independent of the 
metal species. Since the corresponding values for SiC are 
not known, and the buffer layer has a large density of 
states the Fermi l evel, w e simply use the values reported 
for metals in Refs 122122]. 

This "Model 2" gives a stronger doping variation than 
the classical electrostatic model, as show in Fig.|4j:. The 
corresponding resistance is also increased for Model 2, 
as shown in Fig. [4^,, giving excellent agreement with ex- 
periment. Indeed, we consider the degree of quantitative 
agreement between "Model 2" and the experimental data 
to be somewhat fortuitous. Nevertheless, it is striking 
that, using the best approximations available, the mod- 
ulation of local doping by the step can account for the 
observed resistance, while other mechanisms are all far 
too small. 

In principle there could be additional electronic states 
associated with the step that would change the resis- 
tance; but it is not necessary to assume such states in 
order to explain the resistance. Here we assumed the 
vertical surface of the atomic step to be electrically inert. 
This is reasonable since the SiC (0001) surface is electri- 
cally inert, and the extra states are associated with the 
buffer layer. In addition, our results are relatively ro- 
bust against the uncertainties in the bending stiffness 
and the van der Waals binding energy. For example, if 
we assume an unreasonably large Eg — 80meV instead 
of Eb = 40meV, the detachment length for the largest 
step height decreases from 1.56 nm to 1.41 nm, suggest- 
ing a decrease in resistance of only « 10%. A factor of 
2 change in the assumed bending stiffness would have a 
similar effect. 

As seen in Fig.|3]D, the graphene is almost fully depleted 
of carriers in the detached region. This suggests a simple 
model where the graphene is completely undoped over 
the detachment length Our situation then resembles 



4 



the problem of minimum conductivity, often discussed in 
the literature^. Transport in this regime is mediated 
by evanescent modes, but instead of an exponential de- 
cay the graphene bandstructure leads to a unique linear 
("pseudo-diffusive") behavior^, where 71 « f^ d . This 
represents the upper bound for resistance due to doping 
variations, where the doping goes to zero in the detached 
region, and overestimates the calculations of "Model 2" 
by about 50%. Combining this with our previous result 
that id ~ 1.2h s explains the linear scaling of resistance 
with step height (i.e. 1Z s tep oc h s ) observed in experi- 
ment. 

Conclusions: We examined the structural deformation 
of graphene crossing over a substrate step, and the var- 
ious intrinsic mechanisms that may cause electron scat- 



tering at the step. We found that deformation gives only 
a small electron scattering. For graphene on SiC, where 
the substrate induces considerable doping, the dominant 
mechanism is rather the abrupt variation in potential 
and doping due to detachment of the graphene from the 
substrate as it passes over a step. Our result recon- 
cile with the various experimental observations, i.e. that 
TZstep increases with step density®, step heights®, step 
bunching^!, all of them are manifestations of increasing 
id- 
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